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Abstract: In spring 2012 CERN provided two weeks of a short bunch proton beam 
dedicated to the neutrino velocity measurement over a distance of 730 km. The OPERA 
neutrino experiment at the underground Gran Sasso Laboratory used an upgraded setup 
compared to the 2011 measurements, improving the measurement time accuracy. An in- 
dependent timing system based on the Resistive Plate Chambers was exploited providing 
a time accuracy of ~1 ns. Neutrino and anti-neutrino contributions were separated us- 
ing the information provided by the OPERA magnetic spectrometers. The new analysis 
profited from the precision geodesy measurements of the neutrino baseline and of the 
CNGS/LNGS clock synchronization. The neutrino arrival time with respect to the one 
computed assuming the speed of light in vacuum is found to be 5t u = TOF c — TOF v = 
(0.6±0.4 (stat.)±3.0 (syst.)) ns and 8t p = TOF c -TOF p = (1.7±1.4 (stat.)±3.1 (syst.)) ns 
for and P^, respectively. This corresponds to a limit on the muon neutrino velocity with 
respect to the speed of light of -1.8 x 10" 6 < (v u -c)/c< 2.3 x 10" 6 at 90% C.L. This 
new measurement confirms with higher accuracy the revised OPERA result. 
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1 Introduction 

The OPERA neutrino experiment [1] reported in 2011 the observation of an anomalous 
value of the CNGS muon neutrino time-of-flight which could be interpreted in terms of 
superluminal propagation [2]. Continuing carefully to scrutinize the measurement, two 
unaccounted systematic effects significantly affecting the result were found by the Col- 
laboration. A final result taking into account those effects was reported in [3] where no 
significant difference between the speed of light and the neutrino velocity was observed. 

In order to cross-check the revised result and improve the precision of the measure- 
ment, CERN provided a narrow bunch proton beam between the 10th and 24th of May 
2012. OPERA profited of this new dedicated run to improve its timing system and to 
further demonstrate the validity of the procedure used to correct the previously reported 
measurement [3]. In this new run, both the Target Tracker (TT) and Resistive Plate 
Chambers (RPC) were used independently. These detectors are briefly described in Sec- 
tion 3. We report here the results obtained by each sub-detector separately and their 
combined result. Profiting from the information given by the muon spectrometers, results 
are provided separately for muon neutrinos and anti-neutrinos present in the CNGS beam. 
Indeed, ~4% of the Charged Current (CC) interactions in the surrounding rock and inside 
the detector result from interactions [4]. 
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2 The new CNGS narrow— bunch neutrino beam 



In autumn 2011 a first CNGS narrow-bunch wide-spacing neutrino beam was used in order 
to cross-check the neutrino velocity result previously obtained by OPERA with a statistical 
method applied to the data accumulated with the standard CNGS neutrino beam tuned for 
the v T appearance search in z/„ — > v T oscillation [3]. This test proved that it was possible, 
within a short period of two to three weeks of bunched beam, to obtain similar or even 
better accuracy for the neutrino velocity measurement compared to the statistical analysis 
applied to the data accumulated with the standard CNGS beam over years. 

The 2011 narrow-bunch beam had bunches separated by 524 ns yielding a total inten- 
sity of 1.1 x 10 12 protons per cycle. In the 2012 run with narrow-bunch narrow-spacing 
beam (BB) the bunch separation is only 100 ns in order to increase the number of delivered 
neutrinos per time unit. In each CNGS cycle, lasting 13.2 s, a single extraction delivers 
4 batches of 16 proton bunches. Each 1.8 ns long (RMS) bunch contains ~10 n protons 
providing an intensity 6 times higher than during the 2011 BB run. Fig. 1 shows the proton 
waveform of one SPS extraction and Fig. 2 gives a closer view of 4 proton bunches, two 
at the end of the first batch and two at the beginning of the second one, as measured by 
a Beam Current Transformer (BCT) at CERN. The bunch spacing of 100 ns still allows 
to uniquely identify the proton bunch corresponding to a given neutrino interaction in the 
OPERA detector. 

In two weeks run a total of 1.8xl0 17 protons on target (p.o.t.) was delivered. OPERA 
recorded 104 on-time events, 67 involving the TT detectors and 62 involving the RPC 
detectors. 
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Figure 1. Intensity of the 2012 proton beam 
as a function of time for one SPS extraction 
as recorded by the BCT. 



Figure 2. Zoom of Fig. 1 in the region be- 
tween two proton batches. 



Just before this running period CERN installed a monitoring system, called White 
Rabbit [5], at CERN and at LNGS to monitor relevant timing parameters of common 
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systems and also specific parameters of each LNGS experiment. OPERA only used this 
facility to monitor the delay of the external GPS time signal along 8.3 km of the LNGS 
optical fibre and the frequency of the Master Clock (MC) oscillator which compromised 
the measurements performed in 2011. 

3 The OPERA detector and its timing systems 

The OPERA neutrino detector at LNGS is composed of two identical Super Modules, each 
consisting of an instrumented target section with a mass of about 625 tons followed by a 
magnetic muon spectrometer. Each target section is a sequence of walls filled with emulsion 
film/lead plate modules interleaved with pairs of horizontal and vertical 6.7 x 6.7 m 2 planes 
of 256 scintillator strips composing the TT [6]. The TT allows the location of neutrino 
interactions in the target. This detector is also used to measure the arrival time of neutrinos. 
The scintillating strips are read out on both sides through wavelength shifting fibres coupled 
to 64-channel photomultipliers. Each of the two OPERA magnets is instrumented with 22 
planes of RPC's [1]. Each plane can provide the x and y transverse coordinates for crossing 
particles using horizontal and vertical readout strips. 

The DAQ records the timestamp of the earliest TT photomultiplier signal of the event 
and the earliest hit of each RPC module reaching the readout electronics. A time calibration 
of these two sub-detectors allowed converting the measured time in UTC. The complete 
OPERA timing system is extensively described in [3]. Here, only improvements made for 
the new BB run are presented (faster optical/electrical converter, internal Master Clock 
TDC, RPC Trigger Board timing). 

A faster optical/electrical converter was used in order to avoid the effect of time-walk 
discussed in [3] caused by the low light intensity of the lPPmS GPS signal after 8.3 km 
path along the optical fibre from the external to the underground laboratory. After the 
optical/electrical converter processing, the lPPmS signal has a jitter of 3.2 ns (RMS) w.r.t. 
the input signal. The previous measurements of the neutrino velocity were affected by a 
jitter of ±25 ns due to the lPPmS signal tagging with respect to the uncorrelated internal 
frequency of the 20 MHz MC. As a consequence, the timestamp of all hits recorded in 
the same 0.6 s DAQ cycle (the reset of the cycles is sent to all OPERA sensors every 600 
lPPmS pulses) was shifted by the same quantity. In order to remove this jitter an internal 
TDC of the MC FPGA with an accuracy of ~300 ps was used to tag the arrival of the 
external lPPmS with respect to the start of the internal oscillator time-bins, 50 ns wide. 

Because of the 3.2 ns jitter of the lPPmS signal, hits occurring near the end of the 
50 ns time-bins of the internal oscillator are sometimes incorrectly associated with the next 
time-bin. The corresponding event timestamps are thus shifted by 50 ns. The fraction of 
shifted events (outliers) is 5.2%. These outliers are removed by rejecting cases where the 
TDC time (varying from to 50 ns) exceeds 45 ns. Fig. 3 shows the time distribution with 
respect to the beginning of each DAQ cycle for lPPmS pulses coinciding with the DAQ 
reset sent to one TT channel through an injection capacitor. The timing distribution as 
recorded by the DAQ without TDC corrections is shown in blue, the same distribution 
after correction is given in red. The improvement on the time resolution when using TDC 
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Figure 3. Recorded time for lPPmS pulses injected at the level of one TT channel and coinciding 
with the DAQ reset, before (blue) and after (red) TDC correction. 



corrections is clearly seen. The same distribution was also used to accurately estimate the 
time delay between the MC and the reference TT sensor (the time shift between the two 
distributions is arbitrary). 

The RPC timing system is described in [3]. For these new measurements, in order to 
reduce the systematic errors, an accurate determination of all relevant delays was done. The 
starting points of the neutrino interaction time reconstruction are the RPC hits associated 
with the reconstructed muon tracks. The maximum number of hits for tracks crossing one 
complete spectrometer is 22 (44 hits for tracks crossing both spectrometers). The time 
of the signal formation in the RPC gas was obtained from measurements performed on 
a detector prototype. The delays related to the signal propagation in the readout strips 
up to the front end boards and the signal processing in the readout electronics up to 
the timestamp digitization delays were obtained from a dedicated in situ time calibration 
campaign [7]. The lPPmS pulse of the LNGS atomic clock was sent through the same 
readout chain as the one used for the neutrino data. Values are summarized in Tab. 1. 

The OPERA spectrometers are equipped with a new dedicated timing system using the 
purposely developed electronics for triggering the Precision Tracker (PT) drift tubes: seven 
RPC layers in each magnet are instrumented with timing boards [8], each one discriminating 
the positive polarity signals and forming the OR from 16 read-out strips. Each RPC layer 
is served by 14 timing boards (TB), whose digital output signals are sent through 16 m 
long flat cables to the OR Plane Electronics (OPE) boards. Each OPE board forms the 
OR signal of an entire RPC layer which is further used to trigger the PT [9]. An additional 
single-ended positive polarity output is present in the OPE and it was used for the new 
timing system. 

The UTC time of neutrino interactions with muons crossing the OPERA spectrometers 
is reconstructed in two steps: the UTC time coded in the synchronization signal (lPPmS), 
sent by the ESAT LNGS atomic clock, is acquired, while the sub-ms component is computed 
as a time difference between the event hits registered by the TB and the leading edge of 
the lPPmS measured by means of a TDC. 
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Table 1. Parameters used and related systematic errors on the measurement of St. 



Sub-system 


Parameter 


lis 


Common 


R a col i n a ( -4-9fl t^ry~i\ 

Ddseinie ^^zu cm j 


7^1 978 m 7 


l V \ 1 "1 1 1 > 1 / i V ) 
V /1,/lf, l, lilt til 


H\A/ T 1 t~ T1 (TITOT" /~1 

L VV LJ Li ItitiLI Itiy 


±1.0 


Common 


BCT calibration 


583 7 + 1 


Common 


CNGS-OPERA GPS synchronisation 


2.3 ± 1.7 


Common 


UTC delay 


±1.0 


Common 


Total common systematic uncertainty 


±2.5 


TT+RPC (MC) 


LNGS fibre 


41067 ± 1.0 


For RPC-TB 


LNGS fibre 


41319.5 ± 0.2 


TT+RPC (MC) 


DAQ clock transmission 


7046 ± 1.0 


Common 1 1 


rrbA calibration 


24.5 ± 1.0 


/" 1 . . — . . . rpm 

Common 1 1 


Time response 


49.3 ± 0.9 


1 1 Mctnoa 1 


MC simulation 


9.4 ± 3.0 


1 i Mctnoa 1 


TT total intrinsic lst-hit systematic uncertainty 


±3.3 


TT Method 2 


MC simulation 


±0.7 


TT Method 2 


TT Total intrinsic fi systematic uncertainty 


±1.5 


Common RPC 


Signal formation 


24.0 ±2.0 


Common RPC 


Strip propagation delay (4.54 ns/m) 


±0.1 


Common RPC 


Interconnection board 


2.4 ± 0.1 


RPC Method 3 


FEB delay 


±1.0 


RPC Method 3 


Internal mezzanine delay 


±1.0 


RPC Method 3 


Non-uniformities 


±3.9 


RPC (/i) 


RPC Total intrinsic fi systematic uncertainty 


±4.6 


RPC Method 4 


electronic chain delays 


228.0 ±0.8 


RPC Method 4 


TDC ref. frequency 


±0.1 


RPC Method 4 


TDC meas. ref. time 


±0.7 


RPC Method 4 


TDC integral non-linearity 


±0.3 


RPC Method 4 


RPC Total intrinsic TB systematic uncertainty 


±2.3 



The decoding of the lPPmS UTC time is performed by a custom VME module, called 
slave clock, which also reproduces with a negligible time jitter (0.2 ns) the lPPmS signal 
in NIM format. The module output signals, seven for each spectrometer, are discriminated 
and also reshaped in NIM format. 

Signals from RPC's and the lPPmS are converted to differential ECL and sent to a 
CAEN V767 VME TDC module (800 fj,s range, 0.8 ns Least Significant Bit and integral 
non-linearity lower than 0.3 ns, multi-hits feature). The TDC and the slave clock are 
acquired by means of a CAEN V1718 VME bridge. The acquisition system is triggered by 
a 4/14 majority of the RPC layers. 

To obtain a ns-precision over the full range, the TDC was calibrated with a ppm 
precision. For this purpose a thermally stabilized oscillator producing a 5 kHz reference 
signal sent into a TDC channel was used during the whole data taking period. This way, 
the stability of the TDC oscillator can be monitored event by event, performing also the 
calibration with the required precision. 

4 Analysis 

Each event is classified according to its topology [10]: CONTAINED are events with the 
neutrino vertex in the target sections; SPECTRO are events with the neutrino vertex in 
the spectrometer sections; FRONTMUON and SIDEMUON are events with the neutrino 
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vertex outside the detector and with a muon entering from the upstream and lateral sides 
of the detector, respectively. 

Independent analyses were performed for events involving TT and RPC detectors. 
Since TOF c (time of flight) is computed with respect to the origin of the OPERA reference 
frame which is located beneath the most upstream spectrometer magnet, the time mea- 
surements are corrected for the distance of the events along the beam line from this point, 
assuming time propagation at the speed of light. 

Table 2 summarizes the total number of events recorded, as well as the number of 
events used in the four analysis methods described in the next subsections. Methods 1 and 
2 were applied to TT data and Methods 3 and 4 to RPC data. 



Table 2. Events used by the four analysis methods. 



Event type 


Events 


Method 1 


Method 2 


Method 3 


Method 4 


CONTAINED 


18 


17 


13 


13 


11 


FRONTMUON 


37 


29 


27 


27 


24 


SIDEMUON 


17 


3 


3 


9 


8 


SPECTRO 


20 


10 


5 


9 


6 


Others 


12 














Total 


104 


59 


48 


58 


49 



For events where the muon crosses at least one spectrometer its momentum and charge 
are measured. In 3 such events out of 58 CC events used in this analysis, the muon charge 
is positive. One event is CONTAINED (p^ ~12 GeV/c), the other ones are FRONTMUON 
events (p„ ~46 GeV/c and ~ 79 GeV/c). They exhibit a clear muon track crossing both 
spectrometers and thus all drift tube stations could be exploited for charge and momentum 
reconstruction. The overall charge mis-identification probability for these events is below 
1.5%. Considering also the charm background contribution, the number of expected fake 
positive muons from the whole BB sample is < 0.7 for the 3 observed events. If the analysis 
is restricted only to contained events for which both charm contribution and charge mis- 
identification are lower, the expected number of fake positive muons is < 0.07 for 1 observed 
event. 

4.1 Target Tracker data 

All 104 on-time recorded events were analyzed using the TT hit information, applying the 
same selection procedures as described in [3]. 

The selected data sample includes 63 events. Four of them were rejected by the TDC 
cut to remove outliers (as reported in Section 3). 

The neutrino interaction time and the corresponding time of flight TOF u were com- 
puted using two different methods. The first method (Method 1) relies on the information 
of the earliest hit as described in [3]. The second method (Method 2) exploits the time 
information of all TT hits of the 3D muon track. Using the latter method, 48 CC events 
were selected. Method 2 results in a better time resolution (see last column of Tab. 3) 
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Figure 4. Neutrino time distribution of events selected for the TT analysis using the first hit 
(Method 1), for the TT analysis using all muon hits (Method 2), for events selected for the RPC 
analysis using all muon hits (Method 3) and for events selected for the RPC analysis using the TB 
(Method 4). 

due to the larger amount of information associated with the whole track. However, the 
second method misses NC-like events (Neutral Current). In Fig. 4, the upper left (right) 
plot shows 5t = TOF c — TOF u for the first (second) method, respectively. The final result 
is reported in Tab. 3. 

4.2 RPC data using standard DAQ 

For this method (Method 3), only neutrino interactions producing a clear muon track in 
the RPC's were used. For tracks crossing both Super Modules an independent mean time 
was computed using the hits of the first and second Super Module separately. 

The algorithm used to reconstruct the neutrino interaction time was the same as the 
one used for the TT analysis with Method 2. The total number of events used with 
Method 3 is 58 (see Table 2 for details). The corresponding 5t distribution is shown in 
Fig. 4 (Method 3) and the final result is given in Tab. 3. 
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4.3 RPC Timing Board data 

With this new method (Method 4) the extra information provided by the special TDC 
(see Section 3) complements CNGS events acquired by the standard OPERA DAQ. The 
association of the additional data is performed on the basis of the UTC time readout by 
the slave clock and of the VME TDC data. In each event, accidental RPC hits inside 
the multi-hit 800 //s buffer of the VME TDC are rejected considering only those hits in a 
128 ns window around the acquired trigger time. The remaining hits acquired by the VME 
TDC are then corrected by the calibrated delays (for the position along the read-out strip, 
information from the DAQ is again used) and averaged. The fit procedure is similar to the 
one used for Method 3. 

The intrinsic time resolution computed using muons crossing both spectrometers is 
1.1 ns. This value contributes to the width of the distribution in the bottom right plot of 
Fig. 4 (Method 4) , which is a convolution of the proton bunch spread and the uncertainty 
on the interaction point in external events. 

5 Results 

As discussed in the previous section, four different, correlated results were obtained for the 
neutrino velocity measurements using, respectively: 

1. Standard DAQ, TT detectors and earliest hit timestamps, 

2. Standard DAQ and TT hit timestamps of the reconstructed muon track, 

3. Standard DAQ and RPC hit timestamps of the reconstructed muon track, 

4. Timing Board system and RPC detectors using muon hit timestamps. 

Methods 1 and 2 are statistically correlated since they use the same DAQ system 
and sub-detector. They are also correlated with Method 3 using the same DAQ system. 
Method 4 is almost uncorrelated with the other methods because of the different timing 
system used. The residual correlation arises from the systematic error of the common part 
of the timing chain. Tab. 1 reports the list of parameters and associated systematic errors 
of the different measurements. The results obtained with each method are reported in 
Tab. 3, as well as the statistical and systematic errors. 

In order to obtain a single combined result, Methods 2 and 4 were used since they 
have the smallest statistical and systematic errors and - moreover - are almost completely 
uncorrelated thus providing maximum information. Values obtained for Method 1 were 
used for events when Method 2 could not be used because of the absence of a muon track 
(20 events). A standard combination procedure in presence of correlated measurements 
was then used [11]. The final result is: 

5t u = (0.7 ± 0.4 (stat.) ± 1.6 (syst. — uncorr.) ± 2.5 (syst. — corr.)) ns (5-1) 
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Table 3. Indivitual results for v + v for each of the four analysis methods reported in the text. In 
each row, the number of events, St = TOF c — TOF„, the statistical and systematic errors and the 
time resolution on single events are reported (not including the 1.8 ns RMS of the proton bunches). 
All time values are in ns. 



Method 


Events 


5t 


Stat. 


Syst. 


Resol. 


1 


59 


-2.1 


1.1 


4.4 


9.8 


2 


48 


1.2 


1.0 


3.3 


6.5 


3 


58 


-2.5 


1.8 


5.3 


9.5 


4 


49 


0.6 


0.4 


3.6 


1.1 



Summing in quadrature the systematic errors, separately for u and v contributions, 
leads to the results: 

5t v = (0.6 ± 0.4 (stat.) ± 3.0 (syst.)) ns (5.2) 

and: 

5t p = (1.7 ± 1.4 (stat.) ± 3.1 (syst.)) ns (5.3) 

for neutrinos and anti-neutrinos respectively (the slight difference in the systematic errors 
arises from the different contributions of each method to them) . When comparing 5t u and 
5tp, systematic errors largely cancel out. It was assumed that all the NC-like events result 
from v interactions. Since both results are compatible with zero, a limit on the deviation 
from the speed of light was derived (90% C.L.): 

- 1.8 x 10~ 6 < (v v - c)/c < 2.3 x 10~ 6 (5.4) 

and: 

- 1.6 x 10~ 6 < (v p - c)/c < 3.0 x 10~ 6 (5.5) 

for and v^, respectively. It is pointed out that the statistical error in (5.3) was computed 
according to the RMS of the v + v distribution. The above results are in agreement with 
those obtained by the other LNGS experiments participating to this CNGS BB run [12]. 

6 Conclusions 

In May 2012 a two-week dedicated CNGS proton beam was provided to perform a mea- 
surement of the neutrino velocity. The OPERA experiment after improving its timing 
system, has confirmed the result reported in [3], showing no significant deviation of the 
muon neutrino velocity from the speed of light. The present result is 5t = 0.6±0.4 (stat.) ± 
3.0 (syst.) ns, giving the limit -1.8 x 10~ 6 < (v v - c)/c < 2.3 x 10~ 6 at 90% C.L. During 
the same period, 3 muon anti-neutrino interactions were also recorded. The corresponding 
8t value is 1.7 ± 1.4 (stat.) ± 3.1 (syst.) ns compatible with that obtained with neutrino 
events. 
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